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Abstract 

In this work, we consider the issue of pricing exchange options and spread options 
with stochastic interest rates. We provide the closed form solution for the exchange 
option price when interest rate is stochastic. Our result holds when interest rate 
is modeled with a stochastic term structure of general form, which includes Vasicek 
model, CIR term structure, and other well-known term structure models as special 
cases. In particular, we have discussed the possibility of using our closed form solution 
as a control variate in pricing spread options with stochastic interest rate. 
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1 Introduction 



Spread options have become increasingly important. They give the holders the right to call 
or put the spread value of two underlying assets against a predetermined parameter K as 
the strike price. In particular, the spread options reduce to the so-called exchange options 
when the predetermined strike price K is set to zero. Spread options and exchange options 
can be viewed as options to exchange one underlying asset for another with respect to the 
strike price. They are used in many situations. One typical example is that the option 
holder is interested in exchanging one commodity for another commodity. For instance, 
in oil industries, the prices of crude oil and refined oil differ from each other, and both 
prices are fluctuating considerably in response to the weather, regional stabilities of world 
oil production centers, and other human and natural parameters. Oil companies may deal 
with the situations of price fluctuations using the spread options or exchange options. Spread 
and exchange options have been of considerable interests to both practitioners and theoretical 
researchers [0, @, |3|. 

For the exchange option, a closed form solution for its price is available. The valuation 
of exchange option was first studied by Margrabe[|IJ, based on the option pricing theory of 
Black and Scholes[f|], and Merton|5], |J. The derivation of Margrabe is a PDE approach 
However, in Margrabe's derivation, the risk- free rate r is assumed to be a constant, which 
is far from reality. It is of great interest to both practitioners and theoretical researchers 
to investigate whether closed form solution exists when the interest rate is modeled with 
stochastic term structure. 

This paper investigates how stochastic interest rate will affect the exchange option pric- 
ing. The closed form solution for the exchange option's price is given when we assume a 
very general stochastic process for the interest rate, which includes Vasicek model |12|, CIR 
model[13], affine term structure models[l"3] and other interest rate models as special cases. 



To our knowledge, this is the first time to provide the closed form solution for exchange 
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option pricing while stochastic interest rate is taken into account. 



We also argue that to price a European style spread option of general strike price K, one 
may use our closed form solution as a control variate to reduce variance of simulation when 



doing Monte Carlo pricing|ni] for stochastic interest rate. The closed form result presented 



here shall be of interest to both theoretician and practitioners. 

In our discussion below, we assume an exchange economy populated by risk-averse agents 
with increasing preferences, and all economic activities take place in the time interval [0,T]. 
All the possible outcomes of this economy is denoted by a measurable space (Q, T) where Vt 
is the set of all possible states and JF is a sigma algebra of subsets of Q. Information arrival 
in this economy is described by a filtration \Tt] < t < T} with Tt = J 7 , and the agents 
belief is modelled by a probability measure P defined on (D,, J 7 ). 

In following sections, we will first rederive the closed form solution of exchange option 
using the risk-neutral martingale approach. It is shown that this provides identical result 
to the one given by Margrabe using partial differential equation. Then, we discuss, within 
the framework of martingale measure, how to price the exchange option when interest rate 
is stochastic. It is shown that our result is valid for most general term structure, only the 
correlation coefficients between the interest rate and the underlying assets will affect the 
option price. 



2 Exchange Options 

Exchange options can be defined by using two underlying assets or commodities which are 
closely related. This correlation between the two assets or commodities results from demand 
substitution or the potential for transformation. In general, an exchange option has the 
following payoff: 

max{A(5 1 (T) - S 2 (T)), 0} = [A (S 1 (T) - S 2 (T))} + (1) 
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where A = 1 for a call and A = —1 for a put, and Si(T) and S 2 (T) are the underlying asset 
prices at maturity T. 



2.1 Constant interest rate 



Let us first review the case of constant interest rate. The pricing closed form was given by 
Margrabe via PDE method. In the following, we give a short review within the framework 
of risk neutral martingales || [7], [| |10| . Part of the results will be used in the section of 
discussing stochastic interest rate case. 

Suppose that in the physical probability space P), the prices of two underlying 

assets for an exchange option follow the geometric Brownian motions, that is, 

^i^ = II .dt + a i dW i (t), i = l, 2 (2) 

where yU 1; cr 1 , cr 2 , and p l2 = Corr[dWi, dW 2 ] are all constants. The price C(t, Si(t), S 2 (t)) 
of a (European) call exchange option at time t is then given by 

C(t,S l (t),S 2 (t)) = E? \e~^^ ds [S\(T) - S 2 (T)} + ] (3) 

where Q is the corresponding equivalent martingale measure. We should note here that 
these two processes in (|2|) are defined in the physical probability space P), however, 

the general valuation formula in ([3]) is derived in the risk- neutral probability space (Q, J 7 , Q). 
The relationship between the risk-neutral probability space and the physical probability space 
is the standard one, which is described by the Girsanov transformation^]. 

If the risk- free rate r is constant, then the closed form formula for pricing the exchange 
option was first derived by Margrabe (1978) using the partial differential equation approach. 
However, we can also compute the expectation value in the risk neutral space, and the 
exchange option price will follow. It is shown that the approach gives the pricing formula 
identical to the one derived by Margrabe. In this case, we assume that there is no dividend 
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paying during the option's life. Since the price processes of the two underlying assets are 
governed by @, under the equivalent martingale measure Q with constant risk-free rate r, 
we then have 
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where 



Note that we can write 
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where $(•) is the standard Gaussian distribution function. Therefore, by @, (0), and (§), 
we obtain 

E Q[ e -(T-t) [ Sl ( T ) - S 2 {T)}+] 
= S 1 (t)<!>(r]a 1 VT^t - m) - S 2 (t)<S>(r ] a 1 VT^t -m-b). 

To sum up, for constant interest rate r, the price of a call exchange option at time t, is given 
by 

E Q [ e -r(T-t) [5i(r) _ 5a(T)] +-| = S^t)^) - S 2 (t)$(d 2 ) (9) 
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where 

di = log[gi(*)/g2(t)] 1 / 2 2 2 Pl2 a l a 2 )(T - 1), (10) 

d 2 = dt- yjial + al -2 Pl2 a l a 2 ) (T-t), (11) 

and $(■) is the standard Gaussian distribution function. This is consistent with the result 
of Margrabe derived with partial differential equation approach I]. 

As shown above, the exchange option pricing can also be obtained within the framework 
of the risk neutral measure, consistent with the result obtained by Margrabe, which was 
derived with partial differential equation approach. In this case of constant interest rate, 
we wish to note that the interest rate r does not enter the pricing formula explicitly. This 
special feature motivates us to look into the issue of pricing exchange option when interest 
rate is stochastic in a general form. The next subsection discusses this in details within the 
framework of risk-neutral measures. 



2.2 Stochastic interest rate 

In this subsection, we discuss the issue of pricing exchange options when interest rate is 
stochastic. It will be shown below that the option price closed form can be found for most 
general one-factor stochastic interest rate processes ( i.e. one Wiener process ). Our result 
applies to the cases where one describes the interest rate such as Vasicek term structure, 
CIR term structure. These are special cases of our consideration. 

In the following, it is assumed that we are always working in the risk-neutral probability 
space (Q, J 7 , Q). Each process below is referred to this risk-neutral probability measure Q. 
Now assume that the short rate r also follows a Markov diffusion process, that is, 

dr(t) = fi(r(t),t)dt + a(r(t),t)dW (t). 

Here, we do not specify a concrete interest rate model. All we need to assume is that 
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the interest rate is a Markov diffusion process. The interest rate is correlated with the 
two underlying assets of the exchange option being considered. Assume further that the 
correlation matrix of [dWo(t), dWi(t), dW2{t)] is 

1 Poi P02 
Poi 1 P12 
P02 P12 1 

where p 01 , p 02 , and p 12 are constants. Using Cholesky decomposition, we can obtain the 
above correlation structure by setting 

dW (t) = dB (t) 

dW.it) = Poi dB (t) + y/l^p^dB^t) (12) 
^ dW 2 (t) = Po2 dB (t) + ei2z£^ dBl (t) + Ji-fa-iltoz£^dB 2 {t) 

where B (t), Bi(t), and B 2 (t) are three independent standard Brownian motions. For the 
underlying assets, their prices will follow the processes below: 



log^T)] = A,+a ± / dB 1 (s) 
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And hence, by conditional expectation, we can price a call exchange option as 
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Note that given a sample path of {B (s) : t < s < T}, using the results in (|7]) and (|J), we 
obtain the following: 

E? ([5i(T) - S 2 {T)] + \{B (s) :t<s<T}) (17) 



where 



0"! - 0"2Pl2 



'o\ + 0"2 - 2p 12 tTi<T 2 ' 



A 2 - A 



V + a l - ' 2 Pl2< J l< J 2)(T - t) 

b = J{a\ + al-2~ Pl2 a 1 a 2 ){T -t). 



And it is straightforward to check 
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e A l +^a 2 1 (T-t) = 5< i ^ e <TiP i^-5 CT ?Poi( T -*) ; (20) 

Since dB (s) \T t ~ M (0, T — t) , by (|HJ), • ■ ■ , (plD , we then obtain exchange option 
price with stochastic interest rates: Under the conditions (CI) the prices of two underlying 
assets follow the geometric Brownian motions @, and there is no dividend paying during 
the option's life; (C2) If the risk-free rate r is stochastic, then the price of a call exchange 
option at time t, defined by (^), is given by 



V Wx + 2 Pi2°"i°"2)(r - t) 
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where 

log [S 1 (t)/S 2 (t)} + V( a 2Po2 " <7iPoi) + (ffiPoi " ^2P 02 )^ 
di(x) = — 

V (K + ^2 - 2p 12 (7lff 2 ) - (<TiPoi - ^2P 2) 2 ] (T ~ t) 
+ 2 + CT 2 - 2pi2^1^2) - (o-lPoi - (J 2P02) 2 ] ( T - *), 

d 2 (x) = - ^ \A(^i + °2 - 2p 12 ffiff 2 ) - (o-iPoi - 02PO2) 2 ] ( T - *)> 

$(•) is the standard Gaussian distribution function and //, z/) denote a Gaussian density 
function with mean \i and variance v. 

Clearly, if a±p 01 = cr 2 p 2 (that is, the covariance between processes dr(t) and dSi(t) is 
the same as the covariance between processes dr(t) and cLS^i)), then the pricing formulae 
for exchange options are the same for both the stochastic and deterministic term structures. 
Therefore, it is attempting to argue that we could use this solution as a control variate if 
one wants to do Monte-Carlo simulation to price spread option with nonzero strike price K. 



2.3 When underlying assets pay dividends 

The pricing formula for exchange options above when interest rates are stochastic is derived 
with the assumption that the two underlying assets, such as stocks, pay no dividends during 
the options' life. However, in case of the underlying assets also pay constant or known 
dividends, the question will become the general spread options pricing problem. 

Assume the amount of dividends di(ti), • • • , di(t m ) for the first asset and ^2(^1), • • • , 
d2(sk) for the second asset to be paid at the the dates < ti < ■ ■ ■ < t rn < T and < s 1 < 
■ ■ ■ < s k < T respectively are known in advance. Without loss of generality, we can write 
the dividend streams for both assets by 

di(fi),- • • ,di(t n ), i = l,2 

where n is the number of dividends payment dates for both assets. Therefore, at time t, the 



10 



present values of all future dividends will be 



Y,d % {t )e-^^I [m {t^ , = 1,2, 

and the values of all dividends paid after time t and compounded at the risk-free rate till 
the option's maturity date T is given by 



n T 

^^(t,)e^ r(T) %, T ](t,), < = 1,2. 



By the same argument of Heath and Jarrow (1988), we decompose the capital gain processes 
Gi{t) into the asset price process Si(t) and the dividends streams. Assume further that the 
Gi(t) also follows the geometric Brownian motions, that is, 

dGi(t) 



Gi(t) 



Hidt + aidWift), i = 1, 2. 



Then the capital gain Gi(t) may be written by 

n T 



Si(t) + Di(t) 



where 



A(t) = ^^(t,)e^ r(r) % J ,T](t). 
j'=i 

Since the dynamics of the capital gains processes Gi(t), i — 1, 2, under the martingale 
measure Q, is 

dGi(t) 



Gi(t) 



r(t)dt + <TidWi(t), 



and Gi(0) = $(0) and G^T) = S^T) + D^T), % = 1, 2. Therefore, the risk-neural pricing 
formula will be 



if 



e -I t T r{s)ds [s^T) - S 2 (T)Y 



- / ( T r(s)ds 



[Gi{T) — G 2 (T) — (Di(T) — D 2 (T))\ 



In this case, we will have to employ numerical method to value the option price. 
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3 Conclusion 



In this work, we have discussed the issue of pricing exchange options and spread options. 
Closed form for the exchange option price is provided explicitly when the interest rate is 
stochastic. Our result is valid for most general term structure model of one factor, which 
includes Vasicek model, CIR model, and well-known models as special cases. Our result 
indicates that only the correlation coefficients between the interest rate and the underlying 
assets will affect the exchange option price. In one special case, a completely explicit form 
of the option pricing can be obtained. We have argued that it is possible to use this solution 
as a control variate when doing Monte-Carlo simulation to price spread options for nonzero 
strike price K and stochastic interest rate. 
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